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Abstract 

Conventional transport theory is not really applicable to non-equilibrium 
systems which exhibit strong quantum effects. We present two different ap- 
proaches to overcome this problem. Firstly we point out how transport equa- 
tions may be derived that incorporate a nontrivial spectral function as a typi- 
cal quantum effect, and test this approach in a toy model of a strongly interact- 
ing degenerate plasma. Secondly we explore a path to include non-equilibrium 
effects into quantum field theory through momentum mixing transformations 
in Fock space. Although the two approaches are completely orthogonal, they 
lead to the same coherent conclusion. 
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I. INTRODUCTION 



In the past two decades the interest in non-equihbrium quantum systems has grown 
tremendously. It is therefore natural, that more and more examples have been found 
for which the applicability of traditional transport theory is doubtful, i.e., macroscopic 
Boltzmann-equations or Vlasov-Uehling-Uhlenbeck equations fail to describe the system 
adequately. Indeed some of these examples, like dynamical features of the early universe, 
collective effects in fusion plasmas, ultrashort phenomena in semi-conductors, nuclear mat- 
ter in ultrarelativistic heavy-ion collisions and plasma drops composed of quarks and gluons 
have two things in common that require a new paradigm for their theoretical understanding. 

The first of these common factors is the (space and time) density of interactions between 
the components of the system. Some of them are interacting weakly in the sense of having 
a small coupling constant, e.g., as between electrons and photons. However, even for weak 
coupling, a system of many particles experiences synergetic effects like collective participa- 
tion in movement. As a consequence, the systems listed above have excitation spectra which 
differ significantly from those of rarefied gases of free particles. The latter however is a basic 
requirement for the applicability of traditional transport theory. 

The second common factor is the importance of the non-equilibrium aspect on the level 
of the system components: Inhomogeneities occur on spatial and temporal scales that are 
comparable to the intrinsic scales of the physical problem. Particularly interesting examples 
are ultrarelativistic heavy- ion collisions, e.g. between lead nuclei ^osp]^ 160x208 GeV 
laboratory energy. Here, time scales are fm/c and characteristic sizes are a few fm - not 
significantly different from the diameter of the nucleons which compose the initial nuclei. 
Such a significant difference however is another requirement for the description of the system 
in terms of macroscopic equations: Traditional transport theory is correct only to first order 
of gradients in density, temperature etc. 

We may summarize the two common aspects of the systems listed above as 1. the break- 
down of the quasi-particle picture, and 2. the non-separation of space-time scales. The 
necessary theoretical description therefore must account for these two aspects. 

With the present paper we wish to contribute to this new paradigm from two sides: In 
section II we shall point out a way to incorporate nontrivial excitation spectra into transport 
theory (which already contains some non-equilibrium features), and in section HI we will 
clean a path for the inclusion of non-equilibrium effects (gradients in external parameters) 
into quantum field theory. These two completely different approaches lead to the same 
consistent conclusion, summarized in section IV. 

II. FROM STANDARD TO QUANTUM TRANSPORT THEORY 

We perform this step by studying the reverse direction, i.e., starting from a Schwinger- 
Dyson equation for the full two-point function of a fermionic quantum field we derive trans- 
port equations. In contrast to the standard treatment however, we do not perform the usual 
quasi-particle approximation. 

As has been pointed out by various authors, (see ref. for an overview) the description 
of dynamical quantum phenomena in a statistical ensemble necessitates a formalism with a 
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doubled Hilbert space. For our purpose the relevant content of this formalism is, that its two- 
point Green functions are 2x2 matrix- valued. We prefer the technically simpler method of 
thermo field dynamics (TFD) 0, but will keep our derivation sufficiently general to perform 
similar calculations in the Schwinger-Keldysh, or Closed- Time Path (CTP) formalism |^. 

Within this matrix formulation, we consider the Schwinger- Dyson equation for the full 
quark propagator S = So + So(DJ2qS, where 5*0 is the free and S the full two-point Green 
function of the quark field, S is the full self energy and the generalized product © is a matrix 
product (thermal and spinor indices) and an integration (each of the matrices is a function 
of two space coordinates). Throughout this paper we use the convention to write space-time 
and momentum variables also as lower indices, e.g. S^.^ = S(x, y). 

In the CTP formulation as well as in the a = 1 parameterization of TFD 0], the matrix 
elements of S, So and S obey 

oil I c22 _ q12 I q21 sr^ll , ^22 _ ^12 ^21 

•^(0) + '^(0) — '^(0) + '^(0) ) + — —2^ — 2^ . {i.) 

Therefore the four components of the Schwinger-Dyson equation are not independent, the 
matrix equation can be simplified by a linear transformation which one may conveniently 
express as a matrix product [^]. It achieves a physical interpretation only in the TFD 
formalism, see ref. 0]. The transformation matrices B are 

Bw = ( (1 -"'-';) . (2) 

depending on one parameter only. For example, the third term in the Schwinger-Dyson 
equation becomes 

t2( \ Q r^s- Q (Ki ^^-l ( Q T,^ Q something \ , . 

B{n)T3SoQ^Q S{B{n)) =( 5^ © J " 

Here, = diag(l,— 1), T,^'^ are the retarded and advanced full self energy function, and 
gRA gj-Q ^\^Q retarded and advanced full propagator (similarly for Sq) 

= - S'^ , S"" = - . (4) 

The diagonal elements of the transformed equation therefore are retarded and advanced 
Schwinger-Dyson equation. The off-diagonal element is a transport equation. 

We now switch to the mixed (or Wigner) representation of functions depending on two 
space-time coordinates: Exp = Jd'^i^ ~ u) 6xp {}P^{x — yY) T^^y with X = {x + y)l'2., the 
"-sign will be dropped henceforth. The Wigner transform of the convolution S G is a 
nontrivial step: Formally it may be expressed as a gradient expansion 

jd\x-y) exp(iP^(x-y)^) S,, G,, = exp (-iO) txpGxP ■ (5) 

O is a 2nd order differential operator acting on both functions appearing behind it in the form 
of a Poisson bracket OAxpBxp = \ {dxAxpdpBxp — dpAxpdxBxp) ■ We will henceforth 
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use the infinite-order differential operator exp(— iO) = cosO — isinO. Propagator and self 
energy are split into real Dirac matrix valued functions 

= Gxp =F ittAxp ^xp = ReSxp =F iTrr^p • (6) 

Axp is the generalized spectral function of the quantum field. 

Now consider the equations obtained by action of Dirac differential operators (= inverse 
free propagators) on the matrix-transformed Schwinger- Dyson equation ||^. The diagonal 
components are 

Tr [(P^7^ - m) Axp] = cos O Tr [ReExP Axp + Txp Gxp] 

Tr [(P^7^ - m) Gxp] = Tr [1] + cos O Tr [ReS^P Gxp - vr^ Txp Axp] • (7) 

Two important facts about these equations have to be emphasized. First notice that these 
equations do not in general admit a (5-function solution for Axp even in zero order of O. 
In fact, in contrast to other papers ^ we find that there is not such thing as a mass shell 
constraint in quantum transport theory. 

Secondly, the equations do not contain odd powers of the differential operator O. This 
implies, that when truncating the Schwinger-Dyson equation to first order in O (the usual 
order for the approximations leading to kinetic equations), the spectral function Axp may 
still be obtained as the solution of an algebraic equation. 



A. Transport equation 

The off-diagonal component of the transformed Schwinger-Dyson equation reads, after 
acting on it with the inverse free propagator 

with kinetic components = {1 — n) S^"^ + nS'^^ and = (1 — n) + nS^^. Inserting 
the real functions defined before, this leads to a differential equation, which contains all the 
features of our desired quantum transport theory [^,^: 

Tr [(9^7^ + 2 sin O ReS^P + cos O 27rrxp) S^p] = 

2iTr [i sin O p Gxp - cos O Sf p ivr^xp] • (9) 

Note, that here even as well as odd powers of the operator O occur, and the solution in 
zero order O is not trivial. To see this more clearly, we define the generalized covariant 
distribution function Nxp through the equation 

(1 - Nxp) S]^p + Nxp S^^p = . (10) 

It may be proven easily from the Kubo-Martin-Schwinger boundary condition ||^, that with 
this definition of Nxp one reaches the proper equilibrium limit Nxp — > np{E), where np{E) 
is the Fermi-Dirac equilibrium distribution function at temperature T, 

npiE) = (e/^^^-^) + 1)-^ . (11) 
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With this definition, Sxp = 27ri {Nxp — n) Axp, and therefore Nxp is the parameter which 
diagonahzes the the full non-equilibrium matrix-valued propagator through the Bogoliubov 
matrix B from (13) IHI: 



B{Nxp)nSxP {B{Nxp)) 



-1 



G 



XP 



mA 



XP 



Gxp + iTT^xp 



(12) 



Although one may use the generalized distribution function directly in the above eq. @, 
we wish to compare this equation to traditional transport theory. Hence, one more step has 
to be performed, which is the consistent expansion to first order in the operator O. 

As outlined before, we may then use a spectral function which is the solution of an 
algebraic equation. One may argue, that in non-equilibrium states a spectral representation 
of the propagator does not exist in general 0] , but one may still exploit the fact that retarded 
and advanced propagator are by definition analytical functions of the energy parameter in 
the upper or lower complex energy half-plane: 



oo 

S^^^{E, p, X) = ReGxP T ^Axp = j dE' AiE', p, X) 



E-E'±ie' 



(13) 



which is nothing but the Wigner transform of 5"^^ = =F27ri9 (±(xo — yo)) Axy 

Inserting this into the quantum transport equation the yields, correct to first order in 
(see refs. [|],|[ for details): 



Tr 



Axp { - m - ReS^P , A^^p 



iTr 



Axp Nxp^'x'p - {Nxp - 1) ^'xp 



— 1 



, dE , „ 
dr I — sin(rE) Tr 
2tc 



A{X;Po + E, P) 



NxP^'\t + r/2, X- P) - {Nxp - 1) ^'\t + r/2, X; P) 



NJ 



(14) 



In this equation, {■, ■} denotes the Poisson bracket, the index A^ means that the derivatives 
are not acting on Axp- 

The integral over the history of our system constitutes a memory term of this generalized 
transport equation. The memory effect therefore can be attributed to the nonzero spectral 
width of the physical excitations in our system. 

However, even when we send this spectral width parameter to zero in the above equation, 
a new effect remains. Using a simple quasi-particle spectral function 



Axp = Pi^l" -m- ReS^P 



5 



25 {MeV] 
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2.5 5 7.5 10 12.5 

FIG. 1. Time dependent spectral width parameter 74. 
Parameters are 5=0.12, Ti = I MeV, Tf = 200 MeV, m = 10 MeV. 



Thin line: Tt from eq. (R7\), thick line: jt from eq. (19). 



5{Po-E{t,x,P)) 



9((P^ - ReE^p)2 - (M - ReS5^p)2) 



9Po 



(15) 



where i?(t, x, P) is the (space-time and momentum dependent) generalized energy of the 
corresponding "particle" -like state, one obtains after integration over Pq 



dN{x, p, t) ^ dE{t, X, p) dN{t, x, p) dE{t, x, p) dN{t, x, p) 



St 



dt 

N{t,x,p) 



dp 



dx 



dx 



dp 



+ 6St 



N{t,x,p) 



(16) 



The l.h.s. of this equation is the Vlasov (or streaming) part of the standard transport 
equation, it is free of dissipation. 

The first piece on the r.h.s. is the standard "collision integral" of the kinetic equation, 
while 6St is a time-local correction which now contains the three-dimensional Poisson bracket 
of the S^^, S^^ self energies with the function N{t,x,p) 0. We will consider this equation 
in more detail in section III of this paper. 

Before solving the generalized transport equation for a simple example, we would like to 
point out that with our derivation we have indeed taken into account the two topics men- 
tioned in the introductory statements. The breakdown of the quasi-particle approximation 
is accounted for by the nontrivial spectral function, and the non-separation of time-scales is 
accounted for by the Poisson bracket on the r.h.s. of eq. (ffl). 



B. Solution of the generalized transport equation 

In the following, we will concentrate on the memory effects hidden in eq. (0). To 
isolate them from the the gradients in the generalized collision integral 5St [N{t, x, p)] of 
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5 10 15 20 

FIG. 2. Normalized time dependent fermionic distribution function for slow quarks. 
Parameters as in Fig. 1; thin lines: left /npirn^Tf) from the Boltzmann equation (21), right 
Nt/np{m,Tf) from the quantum transport equation (pO|); thick lines: Nj^ /nF{rn,Tf) from the 
generalized kinetic equation (p^). 



eq. ([Tq), we consider a very simple model system, tailored to mimic the thermalization of a 
plasma composed of quarks and gluons. We cannot possibly explain all the relevant physics 



in the present short paper, rather refer to the literature on such systems |]I0[. In view of the 
previous derivations, we are interested in a calculation of relaxation time scales for such a 
plasma, and make some simple approximations: 

1. We assume, that a gas of bosons (gluons) is instantaneously heated to a very high 
temperature. In this gas then eventually quark- ant iquark pairs start to pop up, until 
at the very end a thermal equilibrium in the sense of a degenerate plasma is reached. 

2. We assume, that the self energy function for the quarks is dominated by gluonic 
contributions, and that it does not depend on the energy of the quarks nor on the 
space coordinates. 

3. The gluon background is dominated by external conditions, i.e., we neglect the back- 
reaction of quarks on the gluon distribution. 

4. We neglect the influence of anti-quarks in the spectral function. This restriction is 
removed in an extended version of this application, ref. [p!T|. 



We summarize these assumptions in the following ansatz for the imaginary part of the self 
energy function and for the spectral function of quarks: 

TxP = r, = 7° gT{t) = 7° g (T, e(-t) + Tf e(t)) , (17) 
A{t,E,p) = ^-—^^^^-^—-^. (18) 

TT {E- UJt) + 7t 
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Hence, we approximate the quark spectral function by two time-dependent parameters Ut 
and 7i, which we may interpret as effective energy and effective spectral width. Arguments 
for the validity of this approach are given in ref. |Tl| . 



With the above spectral function the coupled system (|^) reduces to a single nonlinear 
equation for 74, plus the condition ujf = uJq = + . This latter condition is more 
complicated, when the anti-particle piece of the spectral function is taken into account |]lT|. 
The energy parameter is chosen as E = uq, which yields instead of eq. as the Schwinger- 
Dyson equation for the retarded (or advanced) two-point function of the quarks: 

7t = gTi + g{Tf - T,) e(t) (1 - e-27*t j (I9) 

In Fig. 1, the solution of this equations is plotted in comparison to the time dependent 
imaginary part of the self energy function from eq. (|l^. It is obvious, that the solution of 
the nonlinear equation (|T^ approaches the imaginary part of the self energy function with 
a characteristic delay time. In ref. ||Tl| it is discussed how this delay time is calculated from 



the system parameters. 

We now consider three different levels of transport theory for this model, the correspond- 
ing generalized distribution functions relabeled to iVj, and . First of all, due to the 
simplicity of our ansatz we obtain as the full quantum transport equation (|^): 

j^Nt = -2 7i {Nt-nF{m,T{t))) (20) 



with T{t) as defined in eq. ([171). This equation looks surprisingly similar to a kinetic equation 
in relaxation time approach. However, this similarity is superficial: The kinetic equation, 
or Boltzmann equation, derived for our simple model system reads 

jN^ = -2 r, (iVf - n^(m, r(t))) . (21) 

Finally, the generalized transport equation is, correct up to first order in the gradients: 

|iVf = -2r, {N^-npim,Tit))) 

+Ue{t) {g{Tf-Ti)f exp{-2Ttt) 

j^G _ nF{m,Tf)Tf - nF{m,Ti)Ti \ 
' Tf-Ti ) ' ^ ' 

In Fig. 2 we show the numerical solution for Nf, and compare it to the Boltzmann solution 
A^^^ as well as the full quantum transport solution Nt. This comparison of the three methods 
shows, that the full quantum transport equation results in a much slower equilibration 
process than the Boltzmann equation. This result is in agreement with other attempts to 



solve the quantum relaxation problem [p!2| , |l3[ : The quantum system exhibits a memory, it 



behaves in an essentially non-Markovian way. 

In particular, for the physical scenario studied here, the time to reach 1-1/e^ ~ 86 % 



of the equilibrium quark occupation number is almost doubled (14.7 fm/c as compared to 
8.2 fm/c in the Boltzmann case). We furthermore find, that with the generalized transport 
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equation one does at least partially describe the memory effects in a quantum system (the 
characteristic time now is 11.4 fm/c). 

Thus, although we have only used a toy model, it might turn out that quantum effects 
(= memory as described in this contribution) substantially hinder the thermalization of a 
strongly interacting over long time scales. A more thorough discussion of this physical result 
is carried out in ref. I|TT|. 



III. FROM QUANTUM FIELD THEORY TO INHOMOGENEOUS TED 

In the following we describe a system of charged bosons in a non-equilibrium state. For 
simplicity we first consider a space-inhomogeneous situation, but neglect the spectral width 
of the bosons. This view therefore is orthogonal to the previous section, where we ended 
in considering a time-dependent system with translational invariance. While this gave us 
a window to isolate the memory effects of quantum transport theory, we will isolate some 
non-equilibrium quantum features in the present picture. 

As we have pointed out, the description of non-equilibrium systems requires to double 
the Hilbert space. The reason is, that the occupation number, i.e., the interpretation of a 
state as "particle" or "hole" (and hence its temporal boundary condition) may change from 
point to point in order to ensure causality. 

The proper formulation is achieved by expressing each physical (causal) excitation oper- 
ator as the superposition of an operator evolving forward in time and a backward evolving 
operator. Such a superposition is known in the Liouville space, we refer to the literature for 



an introduction |1J|T^. Obviously, there is also a corresponding anti-causal (orthogonal) 
linear combination. Denoting the causal single-particle creation and annihilation operators 
by ttqi, a]^i (the two different charges are distinguished by a lower index / = ±) and the anti- 
causal operators with ~-signs, one may conveniently express this superposition as a matrix 
in three-momentum space [|l^,0: 



aki{t) 



T 



jd'q {Br\t,q,k)y Qty-^^^^ 

Bi{t,q,k)^''-' . (23) 




k is the three-momentum of the modes, therefore in this notation (t) creates a negatively 
charged physical excitation with momentum k, while ak^{t) annihilates a positive charge. 
The operator creates a mode with momentum q and energy Eg = E{q), which 

propagates forward in time. C,qi annihilates a state with the same momentum and energy, 
but propagating backwards in time. These two operators therefore may be combined in a 
Bogoliubov transformation, with a 2 x 2 matrix (compare to the fermion case, eq. (p!2D): 



For simplicity, let us choose a local equilibrium state, i.e., N{t, q, k) is the Fourier transform 
of a space-local Bose-Einstein distribution function 



9 



Ni{t, q, k) = ^ J d'z e-'^'^-'^^ n,it, z, (q + k)/2) 



(25) 



With these creation and annihilation operators we construct two mutually commuting com- 
plex scalar fields (p^, (px see ref. p]J^,p!3[] for details: 



d^k 

7W' 



4_(t) e-"=" + afc+(t)e 



ikx 



al^it) e'*=" + afc+(t) e 



(26) 



Each of these is a representation of the canonical commutation relations, and they may be 

. T 

combined in a statistical doublet 



Before proceeding we would like to emphasize that the above heuristic formulation has led 
us to the creation and annihilation operators known from thermo field dynamics (TFD, f^). 
However, the construction of two commuting representations in non-equilibrium field theory 
is also achieved in other formulations ^ , and their existence has been noticed independently 



by several authors JT^JTS 



A. Effective interaction and diffusion 



In this subsection we will rely more on results obtained in the TFD formalism than in the 
other parts of this paper. However, we believe that the approach described in the following 
is also valid and useful for other methods used in thermal field theory. 

First of all we note, that due to the introduction of the momentum mixing terms the 
equation of motion for the physical creation/annihilation operators is 



d [ aki{t) 



with a "bare Hamiltonian" 



akiit) 



H+Q 



and a second part that vanishes for homogeneous systems: 

T 



Q = Y^ jd^kd^q 

i=± 



-akiit) 



1 1 
1 1 



{-idt + Ek-E^) Ni{t,k,q) 



aqiit) 



(27) 



(29) 



Consequently, by introducing the above linear combinations for the physical creation and 
annihilation operators, the physical fields are no longer free. The time evolution acquires a 
term mixing causal and anti-causal field. 
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Let us emphasize at this point, that one should think of H + Q a.s the Liouville operator 
of our quantum system, which is the generator of the time evolution for the density matrix. 
Its spectral properties then are immediately obvious |I^Jl5| , p!8|JT9[] , e.g., it is not bounded 
from below. 

The mixing part Q of the full Liouvillean vanishes when Ni{t, k, q) is time independent 
and proportional to 6{k — q), which according to ( ]25|) is the case when ni{t,z,p) does 
not depend on the space-time coordinates {t,z). Consequently Q couples the systems to 
gradients in the function ni{t,z,p). 

To see this more explicitly, we henceforth introduce a "trivial" dispersion relation Eg = 



q2 _|_ ^2 g^j^j obtain 



iE,-Eg) N,it,k,q) =-i^^9__ ld'ze-'^'^-''>V.n,it,z,Q) + 0{Vl) 



(30) 



where Q = {q + k)/2. One may argue, that this is an ad-hoc introduction of a coupling 
into the system dynamics. However, this formulation has a deeper foundation: Consider 
for a moment an equilibrium system, because then the Liouville operator is identical to 
eq. (I 



and invariant under the thermal Bogoliubov transformation as ( alak — a\ak 



^t^k — ^t^kj ■ This invariance constitutes of a continuous symplectic symmetry of the 
Liouville space - and such a global symmetry can be made local: We search for those 
operators ^, which diagonalize the Liouville operator locally in space and time. Such a 
requirement connects the parameter ni{t,z,p) of the Bogoliubov transformation with the 
system dynamics. 

Therefore, this process of gauging the symplectic symmetry is entirely equivalent to the 
space-time local Bogoliubov transformation diagonalizing the single-particle propagator that 
was carried out in sect. II. 

On the other hand, it is well known that gauging of a symmetry automatically introduces 
the coupling of the field (p^ to gradients in an external scalar field. Here, this external scalar 
field is the distribution function "field" ni{t,z,p), it is classical and does not possess any 
dynamical feature. To summarize these arguments: The new gradient terms in the time 
evolution occur, because we are implementing causal boundary conditions locally in space 
and time. 

Let us now turn to the generalized transport equation used in the previous section. 
Rewritten for bosons, with an effective energy parameter E that is independent of x it reads 



dni{t,x,k) dE{k) dni{t, X, k) 



dt 



dk 



dx 



St 



ni{t,x,k) 



+ 6St 



ni{t,x,k) 



(31) 



According to our derivation in section II we know that 6St contains a three-dimensional 
Poisson bracket of n and some interaction probability. Obviously this implies that 



d^k 
(2^ 



d-^x 6St 



ni{t,x,k) 







(32) 



and therefore we may write the momentum integral as the three-dimensional divergence of 
a current: 
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(2^ 



ni{t,x,k) 



(33) 



with = d/dx. We use the fact, that the momentum integral over the "standard" coUision 
term St [n{t, x, k)] is zero and therefore obtain for the integrated transport equation 



dtni{t, x) + V,. {ji{t, x) + Jf{t, x)) = 0. 
This is the expression for current conservation, with 



3i{t,x) 



— ni{t,x,k) 



(27r)3 



(34) 



(35) 



as the convective part of this current. Eq. (^) imphes, that in our boson gas a current 
J^i{t,x) arises even for zero "convection" ji{t,x): An inhomogeneous temperature distribu- 
tion gives rise to diffusion. For small gradients of the temperature distribution the diffusion 
current Jf(t,x) will be proportional to \/xni{t,x). The momentum integrated transport 
equation then is nothing but Pick's law. 

More detailed derivations within a non-relativistic field model, as well as exemplaric 
calculations of particle flow and energy flow may be found in ref. 



B. Calculation of the Diffusion Coefficient 

Having outlined the diffusion problem in the previous subsection, we need to generalize 
the momentum mixing to excitations with continuous mass spectrum. This amounts to an 
approximation of the fully interacting quantum fields by generalized free fields as already 
discussed in section II, and will certainly complicate matters very much. 

On the other hand we have outlined above, that only the spatial gradient terms enter 
into the diffusion problem. Hence, for the following we may neglect the time dependence of 
the mixing parameter n. Instead of eq. (|23|) we therefore write 



-akiit) 



JdE jd'qA]'\E,k) {Bi\E,q,k)y f ) e-'^* 

V 'iEql / 



JdE Jd^qAy\E,k) 





# 



T 



1^^' I Bi{E,q,k)^^' 

-^Eql 



where the 2x2 Bogoliubov matrices B have the same form as already given in eq. 
Ni is replaced by 



Ni{E,q,k) = j^^ jd^ze~'<^-^~^^ni{E,z) 



ni{E,z) 



^fS{z){E-t,i(z)) _ ^ 



-1 



(36) 
, but 



(37) 



The ,^-operators have commutation relations 
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6iv6{E-E') 6^{k-k') . 



(38) 



Similar relations hold for the ^ operators, all other commutators vanish, see [0. 

The weight functions Ai{E,k) are positive and have support only for positive energies, 
their normalization is 



dEEAi{E,k) 



dEAi{E,k) = Z, 



kl 



(39) 



The principles of this expansion have been discussed in ref. [0, its generalization to 
non-equilibrium states was introduced in ref. |Q. For equilibrium states the combina- 
tion AsiE^k) = A+{E,k)Q{E) - A-{-E, -k)Q{-E) is the spectral function of the bo- 
son field (px and the limit of free particles with mass m is recovered when AB{E,k) — > 



sign(E) 6{E^ - k^ 



m 



sign(E) 6{E — luI). As we have seen in the previous section. 



the concept of a local spectral function may be applied to non-equilibrium states up to first 
order in the gradient operator O. 

Using these relations we are then able to express the current response of the quantum 
field to the gradients in the distribution function by using the Heisenberg equation of motion 
for time-dependent operators, eq. (|27|), to obtain 



J\t,x) = i drl jit-T,x),Q 



(40) 



to 



where j{x) = i{(j)l.Vx4>x — 4>x^x4>i) is the current operator. 

As has been pointed out in P|,p!3|, the i-th vector component of the Z-charged currents 
generated by the inhomogeneity of the system is 



d^Q Q(^) 



(27r)3 2Z, 



Ql 



dE {Ai{E,Q)y 



dE'E' 



dn^j,iE',Q,x)dAiiE',Q) 



dxU) 



dQU) 



(41) 



In principle, this equation may be used to calculate the diffusion constant for interacting 
boson fields. The result of such a calculation made for pions coupled to nuclear matter 
has led to results which are bigger than the semi-classical transport coefficients p2| , pU[ . by 



factors of 10-100 [p!3|,^. Since the transport coefficients are proportional to the relaxation 
time of the system, this corresponds to the result obtained in section II: Quantum effects 
lead to a slowdown of the relaxation process. 

In order to establish our approach more safely, we now insert a boson spectral function 
with a constant spectral width. However, in contrast to the fermionic case discussed in 
section II, we may then no longer neglect the anti-particles. Hence we use 



AB{E,k) 



2E^B 



TT 



(E2 



El 



7|)' + 4E^l% 



(42) 
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Similar to the fermionic spectral function (|T8D , one obtains this form by summing over four 
simple poles in the complex energy plane, each with the same distance from the real ii^-axis 
and with the same residue. A short calculation then gives the approximate results 



[dE (ME, Q)f = ^ + + 0{^b) 



^Qi 

dn,{E',Q,x)dME',Q) _Q(^) d d 

^ -^^) dQ(^ - TqT ^ ^ ^ ' ^ ^ 

such that an integration by parts gives, up to first order in the gradients of the distribution 
function, 

Ji{x) = - + V. ni{t, x) = -D V. ni{t, x) . (44) 

The diffusion coefficient D therefore diverges when jb ~^ 0. In traditional calculations 
of transport coefficients this diffusion coefficient is obtained as = r {v"^) /3, i.e., as the 
product of relaxation time and square average of the particle velocity . As we have shown 
before, the relaxation time of a non-equilibrium quantum system is infinite if the particles 
do not have a spectral width - and therefore our result has the correct free-particle limit. 

For more detailed considerations however, we would have to define a model for an in- 
teracting system, which e.g. yields the spectral width parameter as function of the system 
characteristics (temperature, chemical potential etc.). This would overstress the goal set 
with the present paper. 



IV. CONCLUSION 

With the present paper we have demonstrated, how one may apply modern techniques 
of "thermal" field theory to non-equilibrium states. This has become important because 
of the quantum character of non-equilibrium systems studied today: A definite breakdown 
of the quasi-particle approximation and the emergence of space-time inhomogeneities on 
scales comparable to those of system components are the reasons for the non-applicability 
of "standard" transport theory. 

We have started with the derivation of a generalized transport equation, which accounts 
for the breakdown of the quasi-particle approximation. In Fig. 2 we have compared the 
solution of this equation with the full Schwinger-Dyson equation as well as with standard 
transport theory. We find, that the nontrivial spectral function one encounters when going 
beyond the quasi-particle picture leads to strong memory effects in transport phenomena. 
In particular, we studied a system which exhibits a characteristic delay time. The memory 
effects are partially taken into account by the generalized transport equation - but it turns 
out, that gradients in the external parameters are equally important to describe the physical 
systems listed in the introduction to this work. 

In the second major part of this work we therefore pointed out, how one may incorpo- 
rate gradient effects into "thermal" field theory. This is achieved by introducing a mixing 
of operators in momentum space. Firstly we have shown that such an approach indeed 
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leads to relaxation currents in inhomogeneous systems. Secondly we have used the mixing 
transformation to achieve a linear response result for transport coefficients. 

This transport coefficient calculation was not pursued numerically. However, we may 
rely on previous numerical estimates of transport coefficients achieved with this method . 
Those numerical estimates give a result very similar to the one obtained in the quantum 
transport model of section II: Quantum systems are relaxing slower than anticipated from 
standard transport theory. 

Hence, the two approaches we have demonstrated in the present paper lead to a co- 
herent result, which is also in agreement with more traditional calculations from statistical 
mechanics. We therefore believe that a clear path exists towards an improved treatment 
of quantum systems in non-equihbrium states, and we express our gratitude to our teacher 
Hiroomi Umezawa for leading us on this path. 
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